feasibility study for hydraulic connection in azarshahr plain coastal aquifer and urmia lake, by geochemical ratios and inverse modeling

introduction azarshahr plain is one of the plains adjacent tothe lake urmia that has been faced, during recent years,withover extraction of groundwater and alsoexposed to salinity enrichment. on the other hand, urmia lake with about 5000 square kilometers is the relatively vastlake with hyper-saline water as a rare biosphere resource in iran and the world. this lake has a great economic, social and bio-spherical importance in iran. it has faced withenormous water depletion during the recent years. many studies has been conductedabout climate change, drought, rivers discharge declining and changing in land use so far that increasing in water utilization in agricultural activities have taken into account as the most important factor in water table reduction of urmia lake.however,there is no study about the hydraulic connection between lake and its adjacent aquifers. perception of spatial changes in groundwater components that is resulted from various inputs is a difficult work, especially if chemical composition or boundary conditions change over time, as it is popular in coastal aquifers. utilization of geochemistry in different water resources connection has done in many studies. geochemists generally seek to findchemical ratios such as /  , / , / +  and / +  andand  isotopes as saltwater intrusion indicators. the present study for the first time investigates the hydraulic connection between urmia lake and azarshahr plain aquifer using geochemical ratios and geochemical modeling. materials and methods   urmia lake is a kind of the widest and shallowest lakes in the world and has hyper saline water.a lot of permanent and seasonal river and streams feed the lake. broad and unwanted changes in water and ecology of this lake have been started about 2 decade before and during this time.thewater volume of the lake decreased toless than 14 billion in 2011from about 32 billionin 1995. azashahr study area is considered as one of the lake twelve adjacent aquifers that is located between 45°,46’ to 45°,50’ longitudinal and 37°,43’ to 37°,52’ latitudinal. its total area is about 457 km  that its plain has an area about 124 km. the highest and lowest height in the study area is 3700m and 1282m, respectively. its average annual precipitation is about 221.2 mm whereas the average annual evaporation is about 1579 mm. and themost important stream in azarshahr plain is azarshahrchai.ithas a southeast- northwest trend and is eliminated to the lake because of wide agricultural usage. on the other hand, the total annual discharge ofthe aquifer is about 90.64 million  that is one of the groundwater depletion and saltwater intrusion factors in the study area. groundwater salinity is increasedin azarshahr plain aquifer because of wide discharge in one side and decrease inurmia lake water level on the other hand. this probably shows that saltwater intrusion to azarshahrplain aquifer is in relation with water level decreasesin urmia lake and consequently reveals the hydraulic connection between them. to investigatethis connection, groundwater sampling and analysis were conductedand hydraulic connection evaluation was done via three steps as below. graphical methods three popular piper, stiff and pie chart graphical method was selectedand the samples were plotted after the analysis and were compared to urmia lake water samples. star diagram also was plotted for the samples and their map was drawn. on the other hand, pie chart that shows the percentage of water major components was calculated and their distribution in study area was plotted and compared to urmia water component percentage. geochemical ratios general status of saline water is  and  dominance with 0.86 molar ratio and of course  exceed over the alkali metals ( and ). also,  is more than  absolutely (i.e. / =4.5-5). on the other hand, continental fresh water is known with its variable composition andits dominant anions are, , and . geochemists try to find alsothe below geochemical ratios generally for saltwater intrusion resource determination: / , / , / +  and / +  and also some isotopes such as boron, oxygen and hydrogen. in the present study,geochemical ratios were calculated and plotted in study area for saltwater intrusion resource finding and also evaluation of hydraulic connection between urmia lake and azarshahr plain aquifer. geochemical inverse modeling hydro-geochemical modeling is an important tool for recognition of geochemical reactions in groundwater. in the modeling process, geochemical reactions during the flow path such as solution, precipitation, mixing and cation exchange are simulated and compared withfinal water composition. solute aggregation in groundwater can be supposed as a reaction of water with surrounding minerals in porous media. hydro-geologists utilize the water chemical data for classification of various mineral balance situations. in this study, ion activity product (iap) and balance constant (k) and also saturation index (si) were used for probable minerals, composition creation and comparison with probable phases in urmia lake water sample. therefore,because of enormous volume of calculations in balance reactions and also finding the balance constant for reactions, phreeqc computer code has been used and probable phases in groundwater mixed samples and urmia lake water sample were created and compared to each other. these were carried out to find if there is any chemical connection between lake and azarshahr aquifer. the mixed samples were created in three paths near the lake. mixing reaction was done between samples number 18 and 19, 20 and 23 and also between 25 and 22 in porous media presentation.  results and discussion geochemicalanalysis results of groundwater samples and urmia lake indicatedthat dominant cation in azarshar aquifer and urmia lake are and , respectively. also groundwater dominant anions are ,  and during the flow direction, whereas dominant anion in lake water is that is distinct characteristic of marine composition. piper diagram revealed that groundwater type in study area are - -  and - -  whereas urmia lake water is in -  type. the nearest sample to urmia lake water composition is sample number7 that is located out of aquifer domain and is affected by anthropogenic activities probably. also stiff diagram of groundwater samples and urmia lake aren’t similar to each other and show the different sources. /  ratio in saltwater intrusion is lower than marine values (i.e. molar ratio